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1 Introduction 



The study of nuclear structure functions has become a fashionable subject. Apart 
from its intrinsic interest, such analysis has a great impact on the interpretation of 
results from heavy ion experiment. At small values of the Bjorken variable x (< 0.01, 
shadowing region), the structure function F2 per nucleon turns out to be smaller in 
nuclei than in a free nucleon |2] . Several explanations to this shadowing have been 
proposed. 

On the one hand, some models use the fact that, in the rest frame of the nucleus 
the incoming photon splits into a qq pair long before reaching the nucleus, and this qq 
pair interacts with it with typical hadronic cross sections, which results in absorption 
in HI HI ini HI IHl IS] • Thus nuclear shadowing is a consequence of multiple scattering 
which in turn is related to diffraction jHl El E] • This relationship will be developed 
in this paper. Equivalently, in a frame in which the nucleus is moving fast, gluon 
recombination due to the overlap of the gluon clouds from different nucleons reduces 
the gluon density in the nucleus OIISI- These studies have received great theoretical 
impulse with the development of semiclassical ideas in QCD and the appearance of non- 
linear equations for evolution in x in this framework (see ^1 ^1 El Ej and references 
therein; also ^H] for a simple geometrical approach in this framework). 

On the other hand, other approaches 1201 1^ do not address the origin of 
shadowing but only its evolution with InQ^: parton densities inside the nucleus are 
parametrized at some scale Qq and then evolved using the DGLAP evolution 
equations. 

The results from different models usually depend on phenomenological assumptions 
and their predictions (notably for small values of x which are of uttermost importance 
to compute particle production at RHIC and LHC) turn out to be very different. 
For example, concerning the Q^-dependence of shadowing, it can be either constant 
in ini ini m IHl in], or die out logarithmically IIHI EUl 121] or behave as a higher- twist 

naiinj. 

In this paper we will use the relation of diffraction to nuclear shadowing which 
arises from Gribov Theory Reggeon Calculus |i24j and the AGK rules |^. In 
this way we obtain a parameter-free description of nuclear structure functions in the 
shadowing region valid for x < 0.01 and < 10 GeV^, using a model for F2 and 
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1211 HZ] • The same strategy has been used in [TUl [TT] , but our extrapolation to 
smaller x or higher W"^ is more reliable than that of fU] due to the model employed for 
the nucleon; besides, our description is valid for small while that of ^T] applies to 
> 4 GeV^. In Section |21 the model will be described. In Section IHl numerical results 
will be presented together with comparisons with experimental data and with other 
models. In Section HI the model will be applied to calculate the multiplicity reduction 
factors |2H1 12n] relevant to compute particle production in heavy ion collisions at RHIC 
and LHC Finally, the last Section will contain our conclusions. 



2 Description of the model 

We assume that nuclei are made of nucleons in the spirit of the Glauber model. In order 
to relate diffraction on nucleons with nuclear shadowing, we will follow the procedure 
explained in jTU] . The 7*-nucleus cross section can be expanded in a multiple scattering 
series containing the contribution from 1, 2,. . . scatterings: 

aA = cT^l^ + cT^? + ---. (1) 

(T^^^ is simply equal to Aaa^adeon- Let us consider now the first correction to the non- 
additivity of cross sections which comes from the second-order rescattering a)^ . In 
Fig. 1 diffractive DIS is shown in both the infinite momentum frame and in the rest 
frame of the nucleon. In Fig. 2 it becomes clear that the square of such contribution 
is equivalent to a double exchange with a cut between the exchanged amphtudes, a 

(2) 

so-called diffractive cut. To compute the first contribution to nuclear shadowing a)^' , 
which comes from these two exchanges, we need its total contribution to the 7*-nucleon 
cross section, which arises from cutting the two-exchange amplitude in all possible ways 
(between the amplitudes and the amplitudes themselves in all possible manners). For 
purely imaginary amplitudes, it can be shown [2H that this total contribution is 
identical to minus the contribution from the diffractive cut. Thus diffractive DIS is 
directly related to the first contribution to nuclear shadowing. The final expression 
reads 



cr 



d%Tl{h) j^^^ dM 



h,f2 rlrrT^ 



dM^dt 



Fj\^{tmin} ■} (2) 



t=0 

with Ta(&) = dzpA{b, z) the nuclear profile function normalized to 1, / d% TAib) = 
1, and the mass of the diffractively produced system. The usual variables for 
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diffractive DIS: Q^, x, and t, or xp = x / (3, (3 = q^^^p , are shown in Fig. 1. 

Coherence effects, i.e. the coherence length of the qq fluctuation of the incoming 
virtual photon, is taken into account through 

FA{tmin) = J d% Mb^/^t~)TA{b), (3) 

with tmin = —rn^N^'p ^^"i the nucleon mass. This function is equal to 1 at x ^ 
and decreases with increasing x due to the loss of coherence for x > Xcrit ~ {^nRa)~^- 

Let us briefly examine 0. Here the real part of the Pomeron amplitude, which is 
small for the value of the intercept which will be used ^7\, A = ap{t = 0) — 1 = 0.2, 
has not been taken into account. Also it has been deduced under the approximation 
R\ ^ Rlf, so the t-dependence of the 7*-nucleon cross section has been neglected. 

For A > 20 a nuclear density in the form of a 3-parameter Fermi distribution with 
parameters taken from will be employed to compute both T^ib) and Q. For 
2 < A < 20 a Gaussian profile function is used jHT] : 

but, in order to take into account the t-dependence for these light nuclei, we make in 
the computation of form factors Q the substitution 

Ra — 'Ra + Rn, Rn = 0.8 fm. (5) 
Finally, for deuteron the double rescattering contribution has the form 



= -2 ™" dt ""^ dM' ""rnucleon 



FD{t), (6) 

i=0 



~oo JM2 . dW^dt 

mm 

where Foit) = e"*, a = 40 GeV-^. 

The lower integration limit in Q and (P) is = 4m^ = 0.08 GeV^, while the 

upper one is taken from the condition: 

^P = ^ I I < a;Pmax =^ M^^^ = Q ^— Ij , (7) 

with xpmax = 0.1; this value was used in [57] motivated by the fact that the model is 
only valid for -C or xp <^ 1, i.e. a large rapidity gap is required. In our case, 
variations of xpmax by a factor 2 do not affect the description of nuclear shadowing at 
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X < 0.01 but the choice xpmax = 0.1 is convenient as it guarantees the disappearance 
of nuclear shadowing at x ~ 0.1 (see below) as in the experimental data. 

and a:pFS(Q^ xp, (3) is provided by the model 



The relation between — 



dM'^dt 



-.1^1" 

A-K^aeni J-oo dxpdt 
da^,^{Q',xp,P) 



dW^dt 

where the usual factorization has been assumed: 



da^.^{x,Q\M\t) _ da^.^{x,Q^M 



2 

Bt 



dM^dt dM^dt 



(9) 

i=0 

with B = 6 GeV^^ (as in [22] , see the discussion there; this value is slightly smaller than 
the experimental values 7.2 ± l.l(stat.)l|];^(syst.) GeV-^ ^1 at (Q^) = 8 GeV^ and 
6.8±0.9(stat.)li-?(syst.) GeV-^ jSl for photoproduction). Note that 7 ' — - 

t=Q 

can be obtained directly from atot- However, the model for diffraction we are using [27] 
has mainly been tested after integration in t (most available data are integrated in t). 
For this reason, we use the integrated expression together with the experimental value 
of B. While this is legitimate at present values of x, it can lead to an underestimation 
of shadowing at very small x, due to the increase of B with energy^. 

The model in [27] is based on the dipole picture of the photon and contains two 
components. The small-distance (S) component corresponds to transverse distances 
r between the q and the q of the dipole such that r < tq, and the large-distance 
(L) component to r > tq, with tq = 0.2 fm. In each component a quasi-eikonal 
iteration is introduced in order to enforce unitarity. Reggeon and Pomeron exchanges 
are allowed. For diffraction, a third component is used, namely a contribution from 
the triple interaction of Reggeons and Pomerons. This model has been designed to 
describe the small x < 10~^, small or moderate < 10 GeV^ region, and it contains 
the basic ingredients which allow its safe extrapolation^ to very small x or high W^. 

Eq. (j2]) corresponds to the case with only two scatterings. Its extension to include 

higher order rescatterings is model-dependent. We will use two models: a Schwimmer 

-'^Nevertheless, the effect is not too large: we have checked that an increase of B from 6 to 7.2 

GeV"^ produces an increase of shadowing for Pb of at most 10 % at a; = 10^^. As estimates indicate 

an increase <50 % in S for the smallest x we have studied, x = 10^^, the increase of shadowing due 

to this effect would be at most ^ 25 % for these values of x. 

^In order to use the model for larger x, 0.01 < a: < 0.1, we have made some modifications in |27| : 

there, in Eq. (26) /3min in the normalization denominators has been set to 0, and in Eq. (25) the 
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unitarization jSo] which is obtained from a summation of fan diagrams with triple 
Pomeron interactions, 



^rnucleony « % ^ _ I) f ^ Q2^T A{h) ' ^ ' 

and an eikonal unitarization, 

= ^rnucleon / d% ^^^^^^f^^ g2) " [-2(^ " l)7^A(&)/(a:, Q')] } , (11) 

where we use the relation cr^*nucieon = ^^q"^" -^2(3^; Q'^) vahd at small x. Here, F2(x, Q^) 
is the nucleon structure function, taken from j^. Both expressions ffTUj) and (fTTjl . 
expanded to the first non-trivial order, reproduce the second rescattering result 
Eikonal unitarization will produce larger shadowing than Schwimmer, as can be ex- 
pected by comparing the second non-trivial order in the expansion of both expressions. 
Finally, 



Cr7* nucleon J M^^^ dW^dt 

as required by consistency with 



FAit-min) (12) 

t=0 



The shadowing in nuclei is usually studied through the ratios of cross sections per 
nucleon for different nuclei, defined as 

R{A/B) = f ^ . (13) 

In the simplest case of the ratio over nucleon (equivalent to proton at small x where 
the valence contribution can be neglected), we get: 

/^-(A/nucleon) ^ / d% , ^ J^j^ q.^^^^,^ , (14) 

i?-^(A/nucleon) = J d'b _ {l - exp [-2{A - l)T^(6)/(x, Q')] } . 

(15) 

To calculate shadowing in photoproduction, x is no longer a relevant kinematical vari- 
able. Instead we use the 7*-nucleon center of mass energy W"^. 

In our framework shadowing can also be studied as a function of the impact pa- 
rameter b: 

R{A/nndeon)'^\b) = ] , (16) 

1 + {A- l)f{x,Q^)TA{b) 



Reggeon-Reggeon contribution has been ignored. These two changes shghtly modify the description 
of diffraction but we have checked that the agreement with experimental data is as good as in the 
original version of the model. 
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R{A/nudeonY'\b) 



1 



{l-exp [-2{A-l)TA{b)f{x,Q^)]}. 

(17) 



2(A-l)T^(6)/(x,Q2) 



Finally, the region of applicability of our model is the same as that of the model 
for diffraction on nucleon |27j, i.e. small x<0.01 and small or moderate Q^<10 GeV^, 




including photoproduction. 

3 Numerical results 

In our model and in [221 work in the small x region and thus no distinction is made 
between protons and neutrons. Although usually joined with straight lines, our results 
are computed at the same (x) and {Q"^) as the experimental data. For the latter, inner 
error bars show statistical errors, and outer error bars correspond to statistical and 
systematical errors added in quadrature. 

In Figs. 3-6 a comparison with experimental data at small x from E665 EZI 
and NMC [HHl IHOl HO] is presented. As expected, eikonal unitarization produces larger 
shadowing than Schwimmer. The agreement with experimental data is quite reasonable 
taking into account that no parameters have been fitted to reproduce the data. Two 
comments are in order: First, for C/D and Ca/D in Fig. 3 which shows the comparison 
with E665 data, shadowing looks overestimated for x ~ 0.01, while in Fig. 5 which 
shows the comparison with NMC data, it looks underestimated. This corresponds to 
the known difference between the results of both experiments for ratios over D, while 
the compatibility is restored jHH] when ratios are computed over C. Second, from Fig. 
6 it becomes clear that the evolution with in the model is too slow at a; ~ 0.01, a 
problem related with the lack of DGLAP evolution in the model j^Z] (see |10| ITU HT| 
for an application of DGLAP evolution to initial conditions). 

In Fig. 7 a comparison of the results of our model with those of others is shown, 
for Q2 = 3 GeV^ (except the results of [H! which are at = 4 GeV^). It can be seen 
that the results of different models agree within 15 % at x ~ 0.01 where experimental 
data exist, while they differ up to a factor 0.6 at x = 10~^. At this x, our results are 
the lowest ones but roughly agree with those of [TH] and with one set of while the 
results from |2I] are the highest ones, and those of [HI 1121 HB] and the second set of [TT] 
lie in between. Let us briefly comment on these models: In fOl 1^ an initial condition 
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is parameterized at some Qq and then evolved using DGLAP; the initial condition is 
fitted from the comparison of the evolved results with experimental data (see 01] for a 
comparison between these two models). is a model which uses a saturating ansatz 
for the total 7*-nucleon cross section in the proton, which is introduced in a Glauber 
expression for its extension to the nuclear case. In ^T] some parameterization of 
hard diffraction at Q^, which as in the present work gives nuclear shadowing through 
Gribov's Reggeon Calculus, is employed; this nuclear shadowing computed at Ql is 
used as initial condition for DGLAP evolution. In j32] a Glauber ansatz provides with 
the initial condition for DGLAP evolution. Finally, in |^ a non-linear equation for 
small X evolution is numerically solved and used in the nuclear case. In view of 
the differences at small x among different models, a measurement of F2 in nuclei with 
~ 10 % precision would be a sensitive test to discriminate among them. Lepton-ion 
colliders jlHI could provide with such data. 

In Fig. 8 our predictions for the ratios D, He, Li, C, Ca, Sn and Pb over nucleon 
for = 0.5, 2 and 5 GeV^ are given for x > 10~*. Let us notice that our model is 
designed for the small x region and that no antishadowing or any other effects relevant 
for a;>0.1 have been introduced. The disappearance of shadowing at a; ~ 0.1 is a 
consequence of both the coherence effects in (jS)) and the vanishing integration domain 
in (j2)), see ((Zj). In Fig. 9 results in photoproduction for the same ratios as in Fig. 8 
are given for W"^ < 10^ GeV^, together with predictions for the evolution of the ratios 
C and Pb over nucleon with impact parameter b. Values as low as 0.3 are reached for 
central Pb/nucleon. This evolution with centrality is very important to compute the 
corresponding evolution of particle production in nuclear collisions, and could also be 
measured in lepton-ion colliders jlE]. 

As a last comment in this Section, let us discuss about the twist structure of the 
model (i.e. its structure in powers of l/Q"^). In the model of [2Z] the unitarity correc- 
tions to the L-component are all of order On the contrary, in the ^-component 
the unitarity corrections are higher-twist (they can be expanded as a sum of terms, 
each one containing an additional factor l/Q"^ as compared to the previous one). The 
fact that diffraction is related to the unitarity corrections allows to study the l/Q"^ 
behavior of shadowing in this model. In order to keep only the leading-twist contri- 
bution (terms oc l/Q"^ in the cross section) we ignore the higher-twist contribution of 
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the S'-component to the diffractive cross section^. The results are given in Fig. 10. 
One can see that neglecting these terms introduces only a small difference. The fact 
that nuclear shadowing corrections are predominantly leading-twist is not unexpected, 
as the diffractive cross section is also leading-twist for the relevant kinematical region 
(indeed, in the model of |2Z] the S'-component diffraction is almost negligible for small 
and/or large M^). This is also seen in the fact that the ratio of diffractive to 
inclusive cross sections does not show any strong Q^-dependence for large [17] . 
Here a comment is in order: in |2ZI a parameter sq is introduced in x and (5 to control 
the limit 0, so that all the equations are written for a; = x + s^jiyV'^ + Q^), 

/5 = /5 + so/(M2 + These terms could mimic higher-twist corrections. In Fig. 10 
we check that the effect of varying this parameter^ from the original sq = 0.79 GeV^ 
to So = 0.2 GeV^ is also very small. So, we can conclude that the contribution from 
higher- twist terms to the shadowing of F2 is small. In contrast, in [llj a large higher- 
twist correction for the shadowing is claimed. The approach in this reference is very 
similar to ours: the authors also compute shadowing from the diffractive cross section, 
but using the HI parametrization [48j. The shadowing obtained in this way for Qq = 4 
GeV^ is then employed as initial condition for DGLAP evolution, taking the shadowing 
for valence quarks from ^Hj- When evolved downwards to = 3 GeV^ a disagree- 
ment is found with experimental data on the ratio Ca over D. This disagreement is 
attributed to higher-twist contributions. 

4 Multiplicity reduction in nucleus-nucleus colli- 
sions 

In this framework it is also possible to study the reduction of multiplicities in nucleus- 
nucleus collisions jlHl I2H1 12^1 • We will denote the rapidity of the produced system in 
the center of mass frame by y* . Shadowing as a function of the rapidity of the produced 
particle can be computed taking into account the general relation with the diffractive 
variables: 

= iw^) ■ <^'' 

■^Concretely, we ignore the 5-component in (17) of [27], and in (20) of |27] we set the exponential 

containing to 1. 

4 And setting c = in Eq. (27) of EH- 
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Then the factor for reduction of multiphcities at fixed impact parameter b is 

R (h^ /^'g RA{s)RB{b-s) 

TAB{b) 

RA{B){b) is given by the r.h.s. of (fT^ multiphed by TA(B)ip) and with f{x,Q'^) substi- 
tuted by F{s,y*) (see below), and 

TAB{b) = J dh TA{s)TB{b- s). (20) 

(|19|) takes into account the summation of Schwimmer's fan-hke diagrams for the pro- 
jectile and target, which are joined by a single Pomeron whose cut gives rise to the 
produced particle (Fig. 11). It follows from AGK cancellation (23] that this is the only 
contribution of this type (more complicated diagrams with lines joining upper and 
lower parts of the diagram, cancel). This provides the justification for the factorized 
expression (fT^ . which is true even if more general rescattering diagrams are taken into 
account. 

The reduction factor as a function of the rapidity of the produced particles F{s, y*), 
can be calculated in several ways. The first one is using (|T2|) . but with the integration 
limits inspired by the parton model for hard processes: for projectile A (target B), 

XAiB) = ^e^'\ (21) 

V -5 

with y* > for the projectile hemisphere and y* < for the target one, and rriT = 



m? + the transverse mass of the produced particle. Substituting in the general 
relation for M^^^., (|7j), we get 



j^4-2(A(S)) 



q2 ( ^Prnax _ A ^ q2 f ^_PrnoWl^^y* _ \ ^ .32) 

\Xa{B) J \ rriT J 



while M^j„ remains fixed and equal to 0.08 GeV^, and = m\. 

On the other hand, we can also compute the reduction factor from the formulas 



F{s,y*) =47r / dy 



Vmax X (l(J 



PPP 



t=0 



'yn.^n crp{s) dydt 

where cTp{s) is the single Pomeron exchange cross section and ^^^^^ the triple Pomeron 
cross section. Using the standard triple Pomeron formula for the latter, we get 



1 da 



PPP 



ap{s) dydt 



= CAexp(Ay), (24) 

t=o 
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with C 



9pp(0)r-ppp{0) 



, QppiS^) the Pomeron-proton couphng and rpppifi) the triple 



4A 



Pomeron couphng, both evaluated at t = 0. In this case, the same integration limits 
used above correspond to: 



In the calculations we have used [C = 0.31 fm^,A = 0.13] taken from jSU] (used in 
[23 1211 )• ^ value rriT = 0.4 GeV is employed by default (in [211122] the nucleon mass 
rriN was used). The sensibility of our results to variations in ttit will be examined. 

In Fig. 12 our results at y* = are presented for AuAu at RHIC energies and for 
PbPb collisions at the LHC, versus impact parameter. Reductions of multiplicities at 
6 = by factors ~ 1/2 for RHIC and ~ 1/3 for LHC are found, with a clear increase 
of the suppression with increasing energy. In Fig. 13 results are presented for AuAu 
at RHIC and PbPb at LHC for different y*. Finally, in Fig. 14 the variation with rriT 
of the results at y* = for AuAu at RHIC is studied. A reduction of the suppression 
with increasing rriT is seen, as expected. Let us make two comments: First, our results 
for the reduction factors are very similar to the ones estimated in fIE\ I29j. It has been 
shown in [211122] that, when these reduction factors are used to correct the results of the 
Dual Parton Model, one obtains a good description of the RHIC data on multiplicities 
and their evolution with centrality. Thus, our results provide a detailed calculation of 
these reduction factors which confirms the estimations in [211 EH] • Second, our results 
are important in studying particle production in heavy ion collisions. In particular, the 
dependence of the reduction factors on ttit gives the variation of shadowing corrections 
with the pt of the produced particle^. 

5 Conclusions 

In this paper, we have used the relation which arises from Reggeon Calculus and the 
AGK rules, between the diffractive cross section measured in DIS on nucleons and the 
^For reduction factors based on other mechanisms, see |51[I52| . 




(25) 



with M^(^(^)) given by ((221), and 




(26) 
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first contribution (i.e. double scattering) to nuclear shadowing. The next contributions 
have been estimated using two different methods for unitarization. In this way we have 
obtained a description of nuclear shadowing, based on the model of [27j for diffraction, 
which agrees with the existing experimental data without any fitted parameter. The 
model is designed for the region of a; < 0.01 and < 10 GeV^, i.e. small x and small 
or moderate Q^. 

The same method has been applied in fUlE]- In CD]; a model for diffraction jSU] 
has been used that takes into account unitarization effects in an effective manner, so 
the extrapolation to smaller x or larger W"^ is not so reliable as in the full unitarization 
program followed in j2Z|; furthermore, the description of diffraction in the model we use 
is substantially better due to the inclusion in the fits of new, more precise experimental 
data. In [llj, a model for diffraction is used in order to obtain an initial condition 
for DGLAP evolution at Qo = 4 GeV^, so their leading- twist description for nuclear 
shadowing is not valid at small Q^. On the contrary, we develop a model valid for 
the full low region which does not correspond to any definite twist but contains 
contributions from all twist orders. Nevertheless, it turns out that, as discussed at 
the end of Section IHl the leading-twist contribution is the dominant one in our model, 
which is in reasonable agreement with the existing experimental data. Precise data 
on the Q^-dependence of nuclear structure functions should disentangle between these 
two possibilities. The existing data from NMC [lOj can be well reproduced within the 
leading-twist DGLAP evolution jTHj with an appropriate set of initial conditions. An 
extension of our results using DGLAP evolution for large values of is thus a natural 
continuation of our work ^l] . 

In this framework we have also obtained the factor for multiplicity reduction in 
nucleus- nucleus collisions. This factor reaches values ~ 1/2 and ~ 1/3 for central AuAu 
and PbPb collisions at RHIC and LHC respectively. It is therefore a very important 
ingredient for the computation of particle production at these energies which should 
be taken into account together with other possible effects. 

Comparison among models shows differences of a factor 0.6 for the ratio of structure 
functions Pb/nucleon at x = 10~^ and = 3 GeV^. These differences have a large 
impact in the computation of particle production in nuclear collisions at energies of 
RHIC and LHC. They should be testable in future lepton-ion colliders flB] . 
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To conclude, the method which we have followed offers a natural link between the 
measurements of nucleon diffractive structure functions and nuclear shadowing, and 
between the latter and the suppression of particle production in nuclear collisions. In 
this way the study of Low x Physics at HERA is linked to that of nuclear structure 
functions at future lepton-ion colliders and with Heavy Ion Physics at RHIC and LHC 

m 
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Figure captions: 



Fig. 1. Diagram showing diffractive DIS with the corresponding kinematical variables 
in the infinite momentum frame (left) and its equivalence in the rest frame of the 
nucleon (right). 

Fig. 2. Diagram showing the equivalence between diffractive DIS and two exchanged 
amplitudes with a cut between the amplitudes. 

Fig. 3. Results of the model using Schwimmer (solid lines) and eikonal (dashed 
lines) unitarization compared with experimental data versus x, for the ratios C/D, 
Ca/D, Pb/D [3^ and Xe/D ^Tj (filled circles correspond to the analysis with hadron 
requirement and open circles to that with electromagnetic cuts, see the experimental 
paper for more details). 

Fig. 4. Results of the model using Schwimmer (open circles) and eikonal (open 
triangles) unitarization compared with experimental data versus A, for the ratios Be/C, 
Al/C, Ca/C, Fe/C, Sn/C and Pb/C ^38^ at two fixed values of x. 

Fig. 5. Id. to Fig. 3 but for the ratios He/D, C/D and Ca/D 

Fig. 6. Results of the model using Schwimmer (solid lines) and eikonal (dashed lines) 
unitarization compared with experimental data versus Q^, for the ratio Sn/C |^ at 
two fixed values of x. 

Fig. 7. Comparison of the results of our model using Schwimmer (solid lines) and 
eikonal (dashed lines) unitarization for the ratio Pb/nucleon with other models, versus 
X at fixed = 3 GeV^. HKM are the results from [21], Sarcevic from Bartels 
from Frankfurt from yj^ (at Q"^ = A GeV^), Armesto from [S] and EKS98 from 

m- 

Fig. 8. Results of the model using Schwimmer (solid lines) and eikonal (dashed lines) 
unitarization for the ratios D/nucleon, He/nucleon, Li/nucleon, C/nucleon, Ca/nucleon, 
Sn/nucleon and Pb/nucleon versus x aX Q"^ = 0.5, 2 and 5 GeV^. 

Fig. 9. Results of the model for = using Schwimmer (solid lines) and eikonal 
(dashed lines) unitarization for the ratios D/nucleon, He/nucleon, Li/nucleon, C/nucleon, 
Ca/nucleon, Sn/nucleon and Pb/nucleon (upper plot), and for different impact param- 
eters h for the ratios C/nucleon (plot in the middle) and Pb/nucleon (lower plot), versus 
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Fig. 10. Results for — 0.5 (upper plot) and 5 (lower plot) GeV^ using Schwimmer 
unitarization for the ratio Pb/nucleon versus x, of the model without modifications 
(solid lines), without the higher- twist contribution in the short-distance component 
(dashed lines), and without the higher-twist contribution in the short-distance compo- 
nent plus some modification in parameters (dotted lines) to check the sensibility of the 
results, see text. 

Fig. 11. Diagram showing the contribution to particle production in the central region 
in AB collisions. 

Fig. 12. Results of the model for the multiplicity reduction factor versus impact 
parameter b y* — 0, for AuAu collisions at ^/s — 19, 130 and 200 GeV per nucleon, 
and for PbPb collisions at ^/s — 5500 GeV per nucleon, in the parton model-like 
reahzation (solid lines) and for [C — 0.31 fm^,A — 0.13] (dashed fines). 

Fig. 13. Id. to Fig. 12 but for AuAu collisions at ^/s = 200 GeV per nucleon and for 
PbPb collisions at ^/s — 5500 GeV per nucleon, for y* = 1, 2 and 3. 

Fig. 14. Results of the model for the multipficity reduction factor versus impact 
parameter b at y* = 0, for AuAu collisions at ^/s = 200 GeV per nucleon, in the 
parton model-like realization (upper plot) and for [C = 0.31 fm^,A = 0.13] (lower 
plot). In each plot, lines from bottom to top correspond to = 0.16, 1, 2, 3, 4 and 
5 GeV2. 
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